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For an isotropic two-dimensional system, when the temperature is lowered toward T = 0, it 
has been found recently that, for the transition from the normal to the superfluid state in the 
paramagnetic limit, the order parameter describing the Fulde-Ferrell-Larkin-Ovchinnikov phase has 
an increasingly complex structure with contributions from an increasing number of wavevectors. 
This cascade of phase transitions is directly linked to the fact that, due to the rotational invariance, 
all the wavevectors directions which can enter the order parameter are degenerate. We study how this 
cascade of phase transitions is modified when one takes into account the anisotropy arising in a real 
solid state compound. For a simple model of anisotropy with elliptical dispersion relation, we find 
surprisingly that the cascade of phase transitions is not modified and the degeneracy with respect 
to the wavevector direction is still present. When we take into account a deviation with respect to 
this elliptical model, which is treated to first order in perturbation, we find that the degeneracy is 
lifted and that, basically, a single wavevector is favored right at the transition. However when one 
enters the superfluid phase, additional wavevectors come in the order parameter and the cascade of 
transitions is still present, though in a modified form. 

PACS numbers : 74.20.Fg, 74.25. Op 



I. INTRODUCTION 

Experimental progress has stimulated current research in the field of Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) 
phases [1,2], with the proposal of several new compounds as displaying these phases. Indeed anomalies in the heavy 
fermion compound CeCoIns have been attributed very recently to FFLO phases [3]. This corresponds to a standard 
3-dimensional situation. However the quasi 2-dimensional case is also of very high interest [4] since this is the 
geometry where one can hope to find compounds with very high critical fields and there have been also recent claims 
of experimental observations [5,6]. In practice this corresponds to compounds made of conducting planes, with very 
small coupling between these planes. Many organic materials and also high T c superconductors fall in this class. In 
such cases for a field strictly parallel to the planes, the orbital currents can not flow perpendicular to the planes. Hence 
they do not come into play for the determination of the critical field, which is only controlled by the paramagnetic 
effect. This is just the dominant one which is under study in the FFLO phases. Naturally in practice the field is not 
exactly parallel to the planes. This brings additional structure due to the appearance of orbital currents and related 
vortex structures [7-11]. Naturally it is important to control first the pure FFLO situation in order to be able to 
disentangle this quite complex situation. 

On the other hand the theoretical situation has been shown recently to be more complex than expected at low 
temperature. This is in contrast with the case of higher temperature where, below the tricritical point, the order 
parameter is believed to be a simple cosine A(r) ~ cos(q.r) at the transition from the normal to the FFLO phase. 
In the isotropic case, corresponding to a circular Fermi line (that is a cylindrical Fermi surface), the direction of the 
q vector is degenerate. However it was first realized [12] by Shimahara that this degeneracy can be used to lower 
the free energy and that the most favorable state contains more than a single cosine in its order parameter when 
the temperature is lowered. Quite recently this investigation has been pushed to lower temperature [13], down to 
zero temperature, and it was discovered that in this regime one has a cascade of phase transitions, with the order 
parameter of each of the FFLO phases being made of the superposition of an increasing number of cosines when the 
temperature is lowered. 

Once this quite peculiar behaviour of the isotropic case is recognized, it is of interest to investigate what is left of 
it in real compounds. Indeed no actual solid state compound is exactly isotropic, although in many cases this may 
be a rather good approximation. In particular the Fermi lines corresponding to the contribution of the C11O2 planes 
in high T c compounds are actually fairly circular (with the center of the circle being at (±7T, ±7r) in the Brillouin 
zone), which is somewhat surprising. It seems clear at first that anisotropy will lift the degeneracy with respect to 
the q vector found in the isotropic case, except for some trivial degeneracies occuring because of remaining discrete 
symmetries. Hence one expects that the transition from the normal state to the FFLO phase will occur for a well 
defined wavevector qo (and its opposite if parity is conserved), leading again to a A(r) ~ cos(qo.r) order parameter. 
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This is indeed what we will find in the general situation. Technically this results from the analysis of the second order 
term in a Ginzburg-Landau expansion. However it is also clear that wavevectors in the vicinity of qo are nearly as 
favored for the order parameter as qo itself. As soon as one goes into the FFLO supcrfluid phase (for example by 
lowering the temperature), they will become a possible choice. Moreover entering the FFLO phase produces a growth 
of the fourth order term in the Ginzburg-Landau expansion. As shown by Larkin and Ovchinnikov, this is this term 
which produces a coupling between the various possible wavevectors. In our specific case, we know from the analysis 
of the isotropic case that it favors the existence of nearby wavevectors and is ultimately responsible for the cascade 
found in this isotropic situation. Hence we may rightly expect to see the cascade appear when we penetrate into 
the FFLO phase, although in a milder form since wavevectors with directions far away from q will clearly be still 
excluded. 

An attractive and simple model to explore more quantitatively this problem is to take an ellipsoidal Fermi line, 
corresponding to an ellipsoidal electronic dispersion relation. On the other hand we will keep a simple BCS type 
attractive interaction, with a constant matrix element which does not vary over the Fermi surface. Clearly this is 
not the most general situation, since we might consider a variation of this coupling on the Fermi surface. This would 
naturally give additional degrees of freedom to the problem, which could combine in various ways to the results we 
will find. But it is reasonable to investigate first the simplest case, where the coupling is constant. 

Surprisingly we find that such an elliptical dispersion relation does not lift the degeneracy found for the isotropic 
case and that the cascade of phase transitions is not modified. Although the reason for this will appear to be rather 
simple, this is a somewhat unexpected result which is of quite important practical interest. Indeed we mentionned 
above that it is not always the case that a circular Fermi line is a bad approximation for a real compound. It is clear 
that the Fermi lines of many more compounds can be closely approximated by an ellipse. Hence in order to have 
a model which displays generic features resulting from the anisotropy of the dispersion relation, we include a small 
additional term to the ellipsoidal dispersion relation which we consider as a perturbation and treat to first order. This 
corresponds to the fairly large class of compounds where the Fermi line is not so far from being elliptical. This leads us 
to the expected situation where, except for discrete symmetries, a single wavevector is favored just at the transition. 
However we will show explicitely that, in agreement with the analysis given above, the cascade of transition appears 
when one enters the superfluide phase. Hence this cascade is not destroyed by anisotropy, as one might have feared, 
but it is merely modified in the way one could expect qualitatively. 

The paper is organized as follows. In the next section we consider the second order term in the Landau-Ginzburg 
expansion, which is the one for which anisotropy is expected to give the most important modification. After handling 
in the first subsection the case of the purely elliptical dispersion relation, the more complex case of a small deviation 
from this ideal situation is treated in the second subsection. Then the fourth order term in the Landau-Ginzburg 
expansion, necessary to obtain the cascade of transitions, is investigated in the third section and it is shown that one 
may simply use the results already obtained in the isotropic case. Finally in the fourth section we study the resulting 
modification of the cascade. 



As indicated in the preceding section we take the simplest model for the electronic interaction responsible for 
pairing, that is a contact potential V(r) = VS(r), since our purpose is to concentrate on the effect of the anisotropy 
in the electronic dispersion relation. This potential is just formally as in the standard BCS theory. Since we expect 
our transitions to be second order, we will use a Ginzburg-Landau approach and expand the free energy difference 
Q = Q s — Q n between the superconducting and the normal state up to fourth order in the order parameter. Hence we 
write Q = CI2 + ^4 where f^2 and ^4 are respectively the second order and the fourth order term. Our main problem 
is to calculate the second order term in the presence of anisotropy. In terms of the Fourier components A q of the 
order parameter A(r): 



II. SECOND ORDER TERM 




(1) 



it reads: 




(2) 



q 



The general expression of ^(q, M; T) is [13] (with slightly changed notation): 
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n 2 (q, p, T) = 1 - G(k)G(k + q) 



(3) 



where the up spin free electron propagator is given by G(k) = (iw n — £k) _1 and G(k) = (— iu> n — ^k) _1 - Here 
u) n = uj n — ip, with w n = ixT(2n + 1) being the Matsubara frequency and p = — m)/2 is half the chemical 
potential difference between up spin and down spins. Moreover £k = — tp is the kinetic energy measured from the 
Fermi surface for p = and e_p is the corresponding Fermi energy. 

It is again convenient to introduce the corresponding expression for q = 0: 

Ao(/x,T) = i-T^G(k)G(k) (4) 

For p = 0, the standard critical temperature T c0 of the BCS phase corresponds to Aj(0,T c0 ) = and more generally 
Ao(p, T sp ) = gives as a function of p the temperature T sp (p) of the spinodal transition line, where the normal state 
becomes absolutely unstable against the transition toward the standard uniform BCS phase. Above the tricritical 
point this line is identical to the standard second order transition toward the BCS phase, but below this tricritical 
point it is located inside the domain for the superfluid phase since it is superseded by the FFLO transition. Hence 
we obtain: 

n 2 (q,p,T)=A Q (p,T) + I(q,p,T) (5) 

with: 

J(q, p, T) = -T^[G(k)G(k + q) - G(k)G(k)] (6) 

n,k 

where the summation over the Matsubara frequencies can be extended to infinity, without any need for the standard 
BCS cut-off, since this summation is fully convergent. As already indicated in Ref. [13], Ao(p,T) is an increasing 
function of p. Since we have fi 2 = at the second order phase transition, we see from Eq.(5) that the highest value 
of p for which this transition will appear corresponds to the minimum of value of p, T). 



A. Elliptical Fermi surface 

We consider first an elliptical dispersion relation: 

£k 2m x 2m y ^ 

for which all the surfaces of equal energies, and in particular the Fermi surface (or precisely the Fermi line, since we 
work in two spatial dimensions), are ellipses. We make no specific assumptions on the effective masses m x and m y , 
so the anisotropy can be very strong. It would seem that this should strongly modify the transition, and in particular 
that the anisotropy will select specific wavevectors q for the space variation of the order parameter. One would expect 
that one of the major axes of the Fermi surface ellipse is preferred. It is easy to see that this is not at all the case. 

Indeed we can convert the summation over k in Eq.(6) into an integral, and then make the change of variables 
k x = \/'>n x K x and k y = y/m y K y , which gives dk x dk y = y/m x m y K x K y , and transforms the electronic dispersion 
relation into 2efc = K^. + Ky = K 2 . If we perform the same change for q, i.e. we set q x = ^m x Q x and q y = ^/rn^Q yi 
we have 2ek+ q = (K + Q) 2 - Hence the integral we find after these changes is formally just the same as for an isotropic 
dispersion relation (except that we have to set the mass m = 1) and we can take up the results of Ref. [13]. This will 
be seen again in the next section when we will calculate this integral in another way. For example this gives at T = 0: 

J(q, p, 0) = ^_JL [ Re ln(l + ^1 - Q*) - In 2] (8) 

where we have set Q — (Q/m)-\Af/2. We obtain that the minimum of /(q, p, 0) is —(In 2/2Tt)^m x m y , which is 
reached for Q = 1, that is e q = p 2 /ep- Therefore the modulus q of the wavevector corresponding to the minimum 
depends on the direction of q, but the minimum of /(q, p,0) itself is the same whatever this direction : we have a 
complete degeneracy with respect to the direction of this wavevector, just as in the isotropic case. Naturally this 
conclusion is also valid at non zero temperature. Accordingly we reach the surprising conclusion that the cascade of 
transitions which has been found for the isotropic situation [13] is not modified when we have an elliptical dispersion 
relation (this will be fully justified below when we will consider the fourth order term). 
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B. Nearly elliptical Fermi surface 



Since we have just seen that an elliptical dispersion relation does not lift the degeneracy with respect to the direction 
of q, we have to go away from this simple dispersion relation and consider a more general case. In order to integrate 
Eq.(6) over k, we will use the fact that the dominant contribution arises actually from the vicinity of the Fermi line. 
This allows to integrate over the variable £ k and to be left with an integration on the curvilinear abscissa Sk along 
the Fermi line Sf- One obtains: 



V G(k)G(k + q) = 7T signer— / ds k |V.e k | 



uj n + § q.Ve k 



(9) 



where the inverse modulus |V.efe| of the gradient of tk is essentially the local density of states. 

Actually we will not consider this general situation, but only the case where the departure from the elliptical 
dispersion relation is small in order to be able to perform explicitely the last integration. We could make a Fourier 
expansion of this small departure along the Fermi line. But in the generic situation the result will be dominated by 
the lowest order relevant Fourier component and we will restrict ourselves to the case where only this component is 
present. Since we will be concerned only by the wavevectors which are near those minimizing our second order term 
f2 2 , this is enough to obtain fairly general conclusions. Therefore this leads us to consider the dispersion relation: 



2r\ 



2m x 2m v m x m v eF 



k 2 k 2 



(10) 



where i/< 1. Depending on the sign of 77, this corresponds to a Fermi line which is, with respect to an ellipse, slightly 
deformed toward a rectangle or toward a rhombus. Lower order Fourier components would only shift the ellipse or 
change its axes. 

It is then simpler to handle this case directly, rather than making use of the above general formula. We make the 
change of variables from (k x , k y ) to (p , ip) defined by: 



k x 
k-u 



lm x p cos ip [1 —p 2 sin 2 ip] 

lm y p sin ip [1 p 2 cos 2 ip] 



(11) 



When this is substituted in Eq.(10), this gives = p 2 /2 to first order in 77. We then make directly this change of 
variable in Eq.(6). The Jacobian is found to be: 



d(k x , ky) 



= y/m x m y p (1 - — p 2 ) 



d{p, ip) 

Since we have tk = p 2 /2, we can again integrate over = — tp. This leads to: 

r 2ir 



If 

V G(k)G(k + q) = 7T ^m x m y (1 - 2rj) signw n — — ^ / 

u ( Z7r ) JO 



dip 



1 



uj n + ^q.Ve k 



(12) 



(13) 



Although the calculation can be carried on at non zero temperature, we will restrict ourselves to the T = case 
in the following. Indeed we want to know how the second order term is modified by a small anisotropy and how the 
cascade of phase transitions is affected. Since this cascade occurs at low temperature, we can compare the effect of 
the small anisotropy, at T = 0, to the effect of a small non zero temperature, evaluated for zero anisotropy. This last 
effect is already known from preceding work [13]. At T = the summation over Matsubara frequencies occuring in 
Eq.(6) can be replaced by an integral, which leads to: 



1 - 2?? 

2tt 



[■2-k 
JO 



dip 
2^ 



In 



1 

Tp 



q-Ve fe 



(14) 



This formula is very similar to the corresponding one obtained by LO [2]. Actually the factor (1 — 2-q) can be 
understood as an effective modification of the density of states. Naturally in Eq.(8) we have to take as it is given 
by the dispersion relation Eq.(10), which leads explicitely to: 
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— q. Vefc = Q [cos <p> + 3r] sin 2(ip + a) sin(<£ + 2a)] 



(15) 



where, together with the change Eq.(ll), we have set q x = ^m x Q cos a and q y — ^/rriy~ Q sm a, with Q = (Q / 'p)\J Tp /2, 
and already made the change of variable <p — > <p + a. 

We are interested in finding the minimum of J(q, p, 0). The situation is here quite analogous to the one found 
without anisotropy. In this last case J(q, p, 0) is given by Eq.(8) and its derivative for its minimum, obtained for 
Q = 1, is discontinuous. It is — oo for Q — > 1_ and is finite and positive for Q — > 1 + , as it can be seen directly. This 
is linked to the fact that the argument of the logarithm in Eq.(14) has a double root for Q = 1. In the presence of 
anisotropy the behaviour of the argument of the logarithm is similar around its minima and we will find again that 
the minimum of J(q, p, 0) corresponds also to the case where this argument has a double root ip m . This implies that 
this argument and its derivative with respect to ip are both zero for ip = p m . To first order in r] the derivative is found 
to be zero for ip m ~ (9r//2) sin(4a). We then request that the argument of the logarithm is zero for <p m . Actually it 
can be seen that, up to first order in 77, this is the same as making this argument zero for p = 0. This leads to the 
conclusion that I(q, p, 0) is minimum for: 

Q = 1 - 3r/sin 2 (2a) (16) 

We are still left with the evaluation of Eq.(14) for this value of Q. It is convenient to take then the derivative with 
respect to i] of the integral in Eq.(14). It can be seen that this derivative is not singular, so we can evaluate it for 
77 = 0. This can be done analytically and gives merely — 3cos(4a). Integrating with respect to i] we finally obtain 
that, to first order in 77, the minimum of I(q, p, 0) is given by (for a fixed direction of q, given by a): 

1 _ 2n 

Min[7(q, p, 0)] = -Jm x m y — — - [In 2 + 37?cos4a] (17) 

This result shows that the degeneracy with respect to the direction of q is lifted as we expected. The minimum of 
7(q, p,0) is found for cos 4a = 1, i.e. a = mr/2, where n is an integer. This corresponds to have the direction of q 
along any axis of our nearly elliptical Fermi line. 



III. FOURTH ORDER TERM 



The fourth order term in the free energy difference is given [2,13] by: 

4 = ^ 4 (qi , q 2 , q 3 , q 4 , T) A qi A* 2 A q3 A; 4 (18) 

with: 

n 4 (qi,qa,q3,q4,M.r) = G(k)G(k + qi)G(k + qi -q 4 )G(k + q 2 ) (19) 

n,k 

and where the q$'s satisfy momentum conservation qi + q 3 = q 2 + q 4 . In the case of an elliptical dispersion relation 
£k = k x /2m x + ky/2m y we can make the same change of variables fcj = y/miKi as for the second order term, and set 
again = ^JralQi. Just as above the integral we find after these changes is formally just the same as for an isotropic 
dispersion relation (with m = 1). Hence, as indicated above, the cascade of phase transitions is completely unaffected 
when the isotropic dispersion relation is transformed in an elliptical one. 

In order to find a modification we have to go away from an elliptical dispersion relation and consider the dispersion 
relation Eq.(10). In order to minimize the free energy for fixed directions of the wavevectors q^ we have first to take 
their lengths in such a way that they minimize the second order term in the free energy, and then to evaluate the 
fourth order term for these lengths. Since in the second term the contributions of the various wavevectors add up, 
as it is seen in Eq.(2), we have just to minimize the contribution of each wavevector. At zero temperature, for a 
given wavevector direction a, the length corresponding to the minimum is given by Eq.(16). However we have to take 
into account that, at non zero temperature, this result is modified. In the isotropic case this modification has been 
obtained in [13] as: 

*--i = |m£ (20) 
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where t = T / p. As we have seen, this result is unchanged for an elliptical dispersion relation provided we replace q 
by Q. We have then to compare the anisotropy modification Eq.(16) to the temperature modification Eq.(20). This 
depends on the angle a for the considered wavevector. However we will be interested only by wavevectors with a 
direction near the one which gives the absolute minimum of the second order term, that is near a = (or a = rar/2). 
The departure from this angle will be typically the angular difference a c between two wavevectors appearing in the 
cascade for the isotropic case, which has been found [13] to satisfy o? c = St ln((27r 3 ) 1 / 2 /t) when t is small. Carrying 
this result into Eq.(16) we see that the temperature dependences of the two modifications are the same, with naturally 
the contribution Eq.(20) being dominant in the limit 77 — > 0, although there is a large numerical coefficient. Anyway 
the temperature dependence of the fourth order term will not be modified by anisotropy. For simplicity we assume 
in the following the small anisotropy situation T] <C 1. In this case we are only left with the dominant modification 
Eq.(20). Accordingly the fourth order term is not modified by anisotropy, and we can make use of the isotropic results 
[13] . ' 



IV. THE CASCADE OF PHASE TRANSITIONS 



In this section we will show that, with the anisotropy that we have considered, the cascade of phase transitions 
found in the isotropic case will also appear, although there are naturally quantitative modifications. Indeed we have 
now a competition between the fourth order term which favors the appearance of the cascade, and the anisotropy 
in the second order term which favors wavevectors direction near the preferred direction. Clearly this makes the 
quantitative study more complicated since we have no longer any symmetry reason to argue that all the wavevectors 
in the expansion Eq.(l) will have equal weight nor any reason to have them with equal angular spacing. Hence a 
detailed quantitative study necessarily involves a good deal of numerical analysis. We will not go in such details since 
the qualitative trend, in which we are actually mostly interested, is on the other hand quite clear. 

The general expression Eq.(18) for the fourth order term in the free energy can be rewritten explicitely [2] as: 

^ = ^£(2-£ Qi , Q J|A Q J 2 |A Q /j(a Qi ^ (21) 

by making use of the fact that all the four involved wavevectors have same length and using momentum conservation. 

The angular dependence of J(otQ ii Q j ) and -/(aQ^Qj) has been already studied [13]. This is the quite singular 
behaviour of J (a) which is responsible for the appearance of the cascade. J(a) is less singular and never plays a 
dominant role. Hence we omit it from now on for simplicity, since it does not modify the qualitative picture. As it has 
been shown in Ref. [13], J (a) is quite large positive at low temperature for a = 0. On the other hand it drops very 
rapidly and becomes negative for some angle a c , which has already been introduced in the preceding paragraph. The 
minimum of J (a) is very near a c . Hence one can view the "interaction" of two wavevectors as being strongly repulsive 
for small angle, and attractive beyond a c with the optimum angle being essentially for a c . This is the physical feature 
which gives rise to the cascade. We will show that it produces the same effect when anisotropy is present. 

For this purpose we will compare the free energy for an order parameter with a single wavevector parallel to the 
direction a — 0, favored by the second order term, to the free energy of an order parameter with three vectors, one 
of them in the direction a — and the two other ones in symmetrical directions with angles ±a c . Actually we omit 
here an additional complication which is again of secondary importance. Indeed it is always more advantageous [13] 
to have wavevectors by pairs with opposite directions. Naturally this is neutral with respect to the second order term, 
which is unchanged when the direction of the wavevector is reversed. On the other hand it is easily seen that this 
is energetically favorable from the effect of the J(a) term, since J(ir) < 0. Moreover the J(a) that we are omitting 
favors also this geometry. As a result the order parameter is a combination of cosines, and for example an order 
parameter with a single plane wave does not occur. However since | J(a)\ is fairly small for a = ir, this brings in 
practice little change in the free energy evaluation, and we omit this complication for simplicity. 

Naturally since we need the effect of the anisotropy to be balanced by the fourth order term, the three wavevector 
phase will appear only when we penetrate in the superfluid domain, but not far from the normal to supcrfluid transition 
since anisotropy is small. We evaluate first the second order term under these conditions. It is still given by Eq.(5). 
The expression of 7(q, p, 0) depends only on the ratio q/p and has already been treated in section II. We have only 
to take into account the fact that A (p,T), given by Eq.(4), depends on the location of (T, p) with respect to the 
transition line. In its vicinity it will vary linearly with the parameters T and p. Since we are at low temperature it 
will depend essentially only on p. Taking the derivative of Eq.(13) with respect to p, in the case q = 0, we obtain: 



Q 2 (q,£,0) = iV 



— 1 + 3 77(1 — cos 4a) 



(22) 



G 



where the evaluation has been made for the optimum value of q for a given angle a. We have taken Eq.(17) into 
account and we have set No = ^Jm x m y ((1 — 2rj)/27r), which is the straightforward generalization of the isotropic 
result. Here p, c is the critical effective field for the normal to (FFLO) superffuid phase at T = 0, for a = 0. 

Calling respectively |A | 2 and |Ai| 2 the weights of the order parameter |A q | 2 for a = and a — a c , we have for 
the free energy of the three wavevector state: 

= n \A \ 2 + 2tt 1 |A 1 | 2 + 1 Jo(|A | 4 + 2|A!| 4 ) + 4Ji|A | 2 |Ai| 2 + 2J 2 |A!| 4 (23) 

where we have set Cio = ^2(0, £t, 0), Hi = fi 2 (a c , /I, 0), Jo = J(0),Ji = J(ct c ) and J 2 = J(2a c ). The minimization 
with respect to |A | 2 and |Ai| 2 leads to: 

VLo + Jo|A | 2 + 4Ji|Ai| 2 = (24) 

2fii + (2J + 4J 2 )|Ai| 2 + 4Ji|A | 2 = (25) 

In practice we can simplify these equations because, as we have seen [13] , we have at low temperature | J\ \ -C Jo and 
I J2I <C Jo- This leads merely to |A | 2 = — do/Jo an d |Ai| 2 = — O1/J0, leading to the free energy: 

n nl + 2nj 

No 2J (2b) 

This is to be compared to the free energy fl/N) = — Og/2Jo for a single wavevector with a = 0, obtained simply by 
setting |Ai| 2 = from the start. This simple result Eq.(26) means merely that the dominant effect is that one gains 
free energy by increasing the number of wavevectors, because one lowers comparatively the effect of the self repulsive 
term Jo in the fourth order term. This is the main origin [13] of the cascade. 

Naturally this occurs only when |Ai| ^ can be satisfied, that is when f2i < 0. This means that the second order 
term is now negative not only for a = 0, but also for a = a c . From Eq.(22) this gives explicitely: 

1 - 4- > 3 77 (1 -cos4a c ) (27) 

For a small departure from an ellipsoidal dispersion relation, i.e. 77 <C 1 this implies that one has to make only a small 
penetration into the superfluid domain to meet the three wavevectors state. Since at low temperature a c is small and 
we have its asymptotic behaviour [13] given above in this regime, we have more explicitely at low temperature for the 
boundary of the domain where the three wavevector state can be found: 

1 - 4- > 2Af 1 a 2 c - 192ryiln((27r 3 ) 1/2 /<) (28) 

fJ-c 

Although the numerical coefficient is large, one has to keep in mind that the right-hand side of this inequality is 
always small since it is necessarily much less than 3 77. One can also see qualitatively how this result is generalized to 
states with more than three wavevectors. For example the state with five wavevectors will appear when the second 
order term becomes negative for the angle 2a c . For small anisotropy, this is again reached for a small penetration into 
the superfluid domain. One can also investigate along the same lines when it is energetically favorable to have two 
wavevector state with angles ±a c /2. With the same simplifications as above, the free energy of this state is found 
to be — 0,y 2 /J where Oi/ 2 = ^2(a c /2, fl, 0). By comparing it to the free energy — Hq/2J for the single wavevector 

with a = we see that the two wavevectors state arises when V2^i/ 2 > f2o> that i s: 

l-j-> 3 t/(1 - cos2a c ) ~ 20.57?a 2 (29) 

Mc v 2 — 1 

Hence this state will appear slightly before the three wavevectors state when one enters into the superfluid domain of 
the phase diagram. 



V. CONCLUSION 



In this paper we have studied the transition from the normal state to the FFLO phases in a two dimensional system 
with anisotropy. More specifically we have considered how anisotropy modifies the cascade of phase transitions, found 



7 



for an isotropic system when the temperature is lowered. We have first considered a simple model with an ellipsoidal 
Fermi line, corresponding to an elliptical dispersion relation. We have found surprisingly that such a modification 
of the dispersion relation does not remove the degeneracy of all wavevectors q on the Fermi line with respect to 
the appearance of a FFLO phase. More generally we have found that, in this case, the cascade of phase transitions 
obtained in the isotropic case is not changed. We have then considered how these results are modified when a small 
term is added to the ellipsoidal dispersion relation, in order to eliminate this degeneracy. We have shown that the 
degeneracy is indeed removed and that, essentially, a single wavevector is selected right at the transition. However 
when one enters in the superfluid phase, a cascade of phase transitions still appears. This shows that, quite generically, 
the cascade found for the isotropic case is modified, but not destroyed, by anisotropy. 
We are grateful to A. Buzdin and C. Mora for stimulating discussions. 

* Laboratoire associe au Centre National de la Recherche Scientifique et aux Universites Paris 6 et Paris 7. 
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